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The stability of MgH^ has been studied up to 20 GPa using density-functional total-energy cal- 
culations. At ambient pressure a-MgHb takes a Ti02-rutile-type structure. a-MgHb is predicted 
to transform into 7-MgH2 at 0.39 GPa. The calculated structural data for a- and 7-MgH2 are in 
very good agreement with experimental values. At equilibrium the energy difference between these 
modifications is very small, and as a result both phases coexist in a certain volume and pressure 
field. Above 3.84 GPa 7-MgH2 transforms into /S-MgHh; consistent with experimental findings. Two 
further transformations have been identified at still higher pressure: i) (3- to <5-MgH2 at 6.73 GPa 
and (ii) 8- to e-MgH 2 at 10.26 GPa. 

PACS numbers: 62.50.+P,61.50.Ks,61.66.Fn 



The utilization of high-pressure technology, which can 
reach the range of giga Pascal (GPa), has made con- 
siderable progress in experimental studies of hydrogen 
storage materials. It has been demonstrated for a num- 
ber of metal— hydrogen systems that application of high- 
pressure is an effective tool to produce vacancies in the 
host metallic matrix, which in turn leads to various novel 
properties |jj . The Mg-based alloys possess many advan- 
tageous functional properties, such as heat resistance, 
vibration absorption, recycling ability, etc. ^). In re- 
cent years, therefore, much attention has been paid to 
investigations on specific material properties of Mg al- 
loys. Magnesium is an attractive material for hydrogen- 
storage applications because of its light weight, low man- 
ufacture cost and high hydrogen-storage capacity (7.66 
wt. %). On the other hand, owing to its high operation 
temperature (pressure plateau of 1 bar at 525 K) and slow 
absorption kinetics, practical applications of magnesium- 
based alloys have been limited. However, it has recently 
been established that improved hydrogen absorption ki- 
netics can be achieved by reduced particle size and/or 
addition of transition metals to magnesium and magne- 
sium hydrides j|. 7-MgH2 often occurs as a by product 
in high-pressure syntheses of technologically important 
metal hydrides like Mg2NiH4. Hence, a complete char- 
acterization of 7~MgH2, in particular, knowledge about 
its stability at high pressure is highly desirable. Present 
theoretical investigation assumes importance as the high- 
pressure x-ray diffraction and neutron diffraction stud- 
ies are unable to identify the exact position of hydrogen 
atoms owing to its very-low scattering cross-section. This 
contribution represents the first theoretical report on the 
high-pressure properties of MgH2 as evident from ab ini- 
tio density functional calculations. 

o>MgH2 crystallizes with Ti02-r-type (r = rutile) 
structure at ambient pressure and low temperature |) . 
At higher temperatures and pressures tetragonal a-MgH 2 
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FIG. 1: Calculated unit cell volume vs. total energy relations 
for MgH2 in actual and possible structural arrangements as 
obtained from VASP. Magnified version of the corresponding 
transition points are shown on right hand side of the figure. 



transforms into orthorhombic 7~MgH2. Recently Bortz 
et al. U solved the crystal structure of 7-MgH 2 (a-PbC>2 
type) on the basis of powder neutron diffraction data col- 
lected at 2 GPa. The a- to 7-MgH 2 transition pressure 
is not yet known. Owing to the difficulties challenged in 
establishing hydrogen position in a metal matrix by x-ray 
diffraction high-pressure information on MgH2 is scarce. 
The aim of the present investigation is to remedy this sit- 
uation by examining MgH2 theoretically at high pressure 
in eleven closely related structural configurations. 

We used the ab initio generalized gradient approx- 
imation (GGA) § to obtain accurate exchange and 
correlation energies for a particular structural arrange- 
ment. The structures are fully relaxed for all the vol- 
umes considered in the present calculations using force as 
well as stress minimization. Experimentally established 
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TABLE I: Optimized structural parameters, bulk modulus (-Bo) and its pressure derivative (B' ) for MgH2 in the different 
structural arrangements considered in the present study. 
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a r = rutile. b To elucidate variations of atomic positions with 
pressure, parameter values are given at equilibrium volume as well 
as at tra n si tio n pressures. c Experimental value Q|. d,e Theoretical 
value |14j, [n5j. f a = anatase. 



structural data were used as input for the calculations 
when available. The full-potential linear muffin-tin or- 
bital (FP-LMTO)Q and the projected augmented wave 
(PAW) implementation of Vienna ab initio simulation 
package (VASP) || were used for the total-energy cal- 
culations to establish phase stability and transition pres- 
sures. Both methods yielded nearly the same result, e.g., 
the transition pressure from a- to 7-MgH2 came out as 
0.385 and 0.387 GPa from FP-LMTO and VASP, respec- 
tively. Similarly almost identical values were obtained 
for ground-state properties like bulk modulus and equi- 
librium volume. Hence, for rest of the calculations we 
used only the VASP code because of its computational 
efficiency. In order to avoid ambiguities regarding the 
free-energy results we have always used the same energy 
cutoff and a similar k-grid density for convergence. In 
all calculations 500 k points in the whole Brillouin zone 
were used for cc-MgH2 and a similar density of k points 
for the other structural arrangements. The PAW pseudo 
potentials (l(J were used for all our VASP calculations. 
A criterion of at least 0.01 meV/atom was placed on the 
self-consistent convergence of the total energy, and the 
calculations reported here used a plane wave cutoff of 
400 eV. 

In addition to the experimentally identified a- and 7- 
modifications of MgH2 , we have carried out calculations 
for several other possible [O types of structural arrange- 
ments for MgH 2 (see Table| and Figs, [j] and | ). The 
calculated total energy vs. volume relation for all these 
alternatives are shown in the Fig. EL In order to get a 
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FIG. 2: The stabilities of various known and hypothetical 
MgHh phases relative to a-MgHh as a function of pressure. 
The transition pressures are marked as a arrow at the corre- 
sponding transition points. 
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clear picture about the structural transition points, in 
Fig. U we have displayed the Gibbs free energy differ- 
ence between the pertinent crystallographic structures of 
MgH2 with reference to a-MgH2 as a function of pres- 
sure. The equilibrium volumes, (30.64 and 30.14A 3 /f.u. 
for a- and 7-MgH2, respectively) are within 1 % of the 
experimental values indicating that the theoretical calcu- 
lations are reliable. The calculated positional parameters 
(Table |) are also in excellent agreement with the experi- 
mental data. The calculated values for the bulk modulus 
(Bo, see Table |) vary between 44 to 65 GPa for the 
various structural arrangements. Among the eleven pos- 
sibilities considered the Ag2Te-type leads to the smallest 
Bq value and the AlAu2-type to the highest (Tabic ||) . 

The calculated transition pressure for the a- to 7- 
MgH 2 conversion is 0.387 GPa (Fig. ||) and as the free 
energy of the two modification is nearly the same at the 
equilibrium volume (Fig. Q), it is only natural that these 
phases coexist in a certain volume range |j . Structurally 
a- and 7-MgH 2 are also closely related, both comprising 
Mg in an octahedral coordination of 6 H which in turn 
are linked by edge sharing in one direction and by cor- 
ner sharing in the two other directions. The chains are 
linear in tetragonal a-MgH2 and run along the four-fold 
axis of its TiOvr-type structure, while they are zig-zag 
shaped in 7-MgH2 and run along a two-fold screw axis of 
its orthorhombic a-Pb02-type structure. The octahedra 
in 7-MgH 2 are strongly distorted. The pressure induced 
a-to-7 transition involves reconstructive (viz. bonds are 
broken and re-established) rearrangements of the cation 
and anion sublattices. The occurrence of a similar phase 
transition in Pb02 suggests that thermal activation or 
enhanced shear is of importance pl|. 

The subsequent phase transition from 7- to /3-MgH2 



FIG. 4: The electronic total and partial density of states for 
a- and 7-MgH2. 

occur at 3.84 GPa. Bortz et al. |(| found no evidence 
for the formation of such a (3 modification up to 2 GPa, 
whereas Bastide et al. Q found that at higher pressure 
(4 GPa) and temperature (923 K) there occurs a new 
phase (viz. (3 in accordance with our findings). How- 
ever, the H positions are not yet determined experimen- 
tally and hence our optimized structural parameters for 
/3-MgH2 should be of particular interest. The volume 
shrinkage at the transition point is 1.45 A 3 /f.u. On 
further increase of the pressure to around 6.7 GPa, (3- 
MgH2 is predicted to transform into (5-MgH2 (orthorhom- 
bic Pbc2\) with a volume shrinkage of 1.1 A 3 /f.u. In 
the pressure range 6.7—10.2 GPa the structural arrange- 
ments /?-MgH2, AuSn2 and <5-MgH2 are seen to lie within 
a narrow energy range of some 10 meV (Fig. |]J). This 
closeness in energy suggests that the relative appearance 
of these modification will be quite sensitive to, and eas- 
ily affected by, other external factors like temperature 
and remnant lattice stresses. A transformation to e- 
MgH 2 (AlAu2-type structure also called cotunnite-type) 
is clearly evident from Fig. ^ with a volume change of 
1.7 A 3 /f.u. at the transition point (see also Figs. || and 
^). It is interesting to note that similar structural tran- 
sition sequences are reported for transition-metal oxides 
like Hf02 |pT| . Recently, the cotunnite-type structure 
of Ti02 (synthesized at pressure above 60 GPa and at 
high temperatures) has been shown |l^] to exhibit an 
extremely high bulk modulus (431 GPa) and hardness 
(38 GPa). Our present study predicts that, it is possi- 
ble to stabilize MgH2 in the AlAu2-type structure as a 
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soft material (Bq = 65 GPa) above 10.26 GPa. Hence, 
the present prediction of stabilization of a contunnitc- 
type structure for MgH 2 at high pressure should be of 
considerable interest. 

The calculated unit-cell volume difference between the 
involved phases at their equilibrium volume relative to 
the a-phase is 0.497 (for 7), 2.94 (for /?), 3.60 (for 5), 
and 5.96 (for e) A 3 /f.u., which is approximately 1.6, 9.5, 
11.8, and 19.5%, respectively, smaller than the equilib- 
rium volume of the a-phase. The energy difference be- 
tween a-phase and the 7, 0, S, and e phases at their 
equilibrium volume is 0.81, 43, 66, and 197 meV/f.u., 
respectively, which is considerable smaller than for simi- 
lar types of phase transitions in transition metal dioxides 
lim . The known stabilization of the high-pressure phase 
of TiC>2 by using high-pressure— high-temperature syn- 
thesis |l2| indicates that it may be possible to stabilize 
the high-pressure phases of MgH2 due to their small en- 
ergy differences and (partly) reconstructive phase tran- 
sitions. Such stabilization of MgH 2 high-pressure phases 
at room temperature would reduce the volume consider- 
ably, and in the extreme case the expected volume re- 
duction would be ca. 19.5% for e-MgH2 compared with 
a-MgH2 . This would imply an increased volumetric stor- 
age capacity of hydrogen for MgH 2 . Of interest would be 
to explore the possibility of stabilizing the high-pressure 
phases by chemical means. At the same time this may 
possibly open up for improved kinetics with respect to 
reversible hydrogen absorption/desorption. Hence, the 
observation of high-pressure phases in MgH2 may have 
technological significance. 

The calculated energy gap (from the top of the va- 
lence band (VB) to bottom of the conduction band) is 
4.2 and 4.3 eV for a- and 7-MgH2, respectively. Accord- 
ing to Fig. I the width of VB is 7.5 and 8.3 eV in a- 
and 7~MgH2, respectively. The increased width of VB 



in 7-MgH 2 is due to the reduction in the Mg— H sepa- 
ration (1.93 A in a-MgH 2 vs. 1.91 A in 7-MgH 2 ). An 
experimental UV absorption study [jlJJ of a-MgH2 gave 
an absorption edge at 5.16 eV, whereas earlier theoreti- 
cal electronic structure investigation |l4|, [l5| reported a 
band gap around 3.4 eV. Hence, our calculated energy 
gap value is in better agreement with the experimen- 
tally measured value than the earlier values. In fact, the 
present underestimation of the band gap by about 17 % is 
typical for the accuracy obtained by first principle meth- 
ods for semiconductors and insulators. Such distinctions 
probably originates from the use of GGA-based exchange 
correlation functionals. 

In conclusion, the calculated structural parameters for 
a- and 7~MgH2 are in excellent agreement with the ex- 
perimental findings. We found that the ground state of 
a-MgH2 becomes unstable at higher pressure and the cal- 
culated transition pressures for a-to-7- and 7-to-/?-MgH2 
are found to be in very good agreement with the exper- 
imental findings. We predict that further compression 
of /3-MgH2 will lead to a phase transition to <5-MgH2 
(orthorhombic; Pbdl\) at 6.73 GPa and that the AIAU2- 
type structure of a hypothetical e modification stabilizes 
above 10.26 GPa. Our total energy study suggested that 
the formation of the high-pressure e-phase could also be 
possible at the ambient pressure/temperature by choos- 
ing an appropriate preparation method. This may reduce 
the volume requirement by ca. 19.5% compared with a- 
MgH2. In agreement with experiments, a-MgH2 turns 
out to be an insulator, and the theoretical treatment 
shows that the high-pressure modifications also should 
exhibit insulating behavior. 
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